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Abstract 

A search for the production and non-standard decay of a Higgs boson, h, into four taus 
through intermediate pseudoscalars, a, is conducted on 683 pb -1 of data collected by the 
ALEPH experiment at centre-of-mass energies from 183 to 209 GeV. No excess of events 
above background is observed, and exclusion limits are placed on the combined production 
cross section times branching ratio, £ 2 = CT ^ 6 ( e + e -^>zh) x — aa ) x -^( a ~~ r+r ~) 2 - F° r 
m h < 107 GeV/c 2 and 4 < m a < 10 GeV/c 2 , £ 2 > 1 is excluded at the 95% confidence 
level. 
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1 Introduction 



Searches conducted at LEP2 have excluded the standard model (SM) Higgs boson 
decaying into bb or t + t~ for masses below 114.4 GeV/c 2 pQ. The LEP experiments 
observed a small excess in the bb final state for a Higgs boson mass around 100 GeV/c 2 , 
which is consistent with background fluctuations or SM-like production with a reduced 
branching ratio into bb [2J [3]. This excess, the mild tension with electroweak precision 
tests [I], and the fine-tuning needed in the minimal supersymmetric standard model 
(MSSM) have prompted the consideration of models with exotic Higgs boson decays, such 
as those of the next-to-minimal supersymmetric standard model (NMSSM) [51 [6] as well 
as more general frameworks [HE]. In these models, new decay channels can dominate over 
h — > bb and render the Higgs boson "invisible" for conventional searches. In particular, a 
Higgs boson decaying into two light pseudoscalars is well motivated by these models and 
results in a four-body final state as the pseudoscalars decay into light fermions. Scenarios 
with a dominant h — > 2a — > 4b decay are constrained by the LEP experiments for mass 
up to 110 GeV/c 2 [21 |H1 [TO]- A general search for h — > aa with a — > gg, cc, r + r _ was 
performed by the OPAL Collaboration [IT] , but the analysis was restricted to a Higgs 
boson mass, mh, in the range 45-86 GeV/c 2 . A search for h 2a ->■ 2/i2r was recently 
reported by the DO Collaboration, which resulted in limits that are a factor of 1-4 larger 
than the expected production cross section assuming branching ratios of the pseudoscalar 
as predicted by the NMSSM [12]. In this paper, a search using ALEPH data is presented 
for h — > 2a — > At up to m h « 110 GeV/c 2 . 

The pseudoscalar a may arise from a two Higgs doublet model, as in the MSSM, or 
it can include a component from an additional singlet field as in the NMSSM. These 
possibilities differ in their details and relations between model parameters. The present 
search is performed in a model-independent manner and simply adopts the two main 
characteristics of the pseudoscalar: the coupling to a Higgs boson resulting in h — > aa 
decay and the coupling to SM fermions proportional to their Yukawa couplings. The 
present analysis concentrates on the region 2m T < m a < 2mt,, where the a — > t + t~ decay 
mode is expected to be substantial. The Higgs boson production mode considered here is 
the Higgsstrahlung process, shown in Fig. [ljwith Z — > e + e~, /i + yU~, vv. 




Figure 1: Higgs boson production and decay modes considered in this analysis. 

This paper is organized as follows. The ALEPH detector is described in Section [2J 
Simulation of background and signal processes are described in Section [3J The event 
selection criteria are described in Section HI Systematic uncertainties are discussed in 
Section |5j Finally, results and conclusions are presented in Sections [6] and [71 respectively. 
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2 The ALEPH Detector 



A detailed description of the ALEPH detector can be found in Ref. [13] and of its 
performance in Ref. [H]. Charged particles are detected in the central part, which 
consists of a precision silicon vertex detector, a cylindrical drift chamber and a large 
time projection chamber (TPC), together measuring up to 31 space points along the 
charged particle trajectories. A 1.5 T axial magnetic field is provided by a superconducting 
solenoidal coil. Charged particle transverse momenta are reconstructed with a 1/pr 
resolution of (6 x 10 -4 @5 x 10~ 3 /p T ) (GeV/c) -1 . The charged particle tracks used 
in the present analysis (and simply called tracks) are reconstructed with at least four hits 
in the TPC, and originate from within a cylinder of length 20 cm and radius 2 cm coaxial 
with the beam and centred at the nominal collision point. 

Electrons and photons are identified by the characteristic longitudinal and transverse 
development of the associated showers in the electromagnetic calorimeter, a 22-radiation- 
length thick sandwich of lead planes and proportional wire chambers with fine readout 
segmentation. A relative energy resolution of 0.18/ a/E (E in GeV) is achieved for isolated 
electrons and photons. 

Muons are identified by their characteristic penetration pattern in the hadron 
calorimeter, a 1.2 m thick yoke interleaved with 23 layers of streamer tubes, together with 
two surrounding double-layers of muon chambers. In association with the electromagnetic 
calorimeter, the hadron calorimeter also provides a measurement of the hadronic energy 
with a relative resolution of 0.85/y/E (E in GeV). 

The total visible energy is measured with an energy-flow reconstruction algorithm 
which combines all the above measurements [14J. The relative resolution on the total 
visible energy is 0.60/ \f~E for high multiplicity final states. In addition to the total 
visible-energy measurement, the energy-flow reconstruction algorithm also provides a list 
of reconstructed objects, classified as charged particles, photons and neutral hadrons, and 
called energy-flow objects in the following. These energy-flow objects are the basic entities 
used in the present analysis. 

Below polar angles of 12° and down to 34 mrad from the beam axis, the acceptance 
is closed at both ends of the experiment by the luminosity calorimeter (LCAL) and 
a tungsten-silicon calorimeter (SICAL) originally designed for the LEP 1 luminosity 
measurement. 

The average centre-of-mass energies at which the machine operated and the 
corresponding integrated luminosities used in this analysis are presented in Table [TJ 

Table 1: Integrated luminosities collected at the different average centre-of-mass energies. 



^CM(GeV) 


182.65 


188.63 


191.58 


195.52 


199.52 


201.62 


204.86 


206.53 


jCdt (pb- 1 ) 


56.8 


174.2 


28.9 


79.9 


86.3 


41.9 


81.4 


133.2 
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3 Signal and Background Generation 



Both signal and background were generated for all centre-of-mass energies shown in 
Table [1] using the GEANT3-based simulation of ALEPH [15]. Backgrounds were generated 
with a variety of generators listed in Table [2j The 77 initiated and Bhabha samples are 
10-30 times larger than the data sample while others are 300-1000 times larger. The 
contribution to the Wez/ process from electrons emitted close to the beam axis is not 
included in koralw and is generated with pythia. The HZHA03 generator [16] was used 
to generate 3000 signal events (with h — > aa followed by a — > t + t~) for each of the three 
Z decay channels considered and for each combination of Higgs boson and pseudoscalar 
masses in the ranges 70 < < 114 GeV/c 2 and 4 < m a < 12 GeV/c 2 in steps of 
2 GeV/c 2 . 

Table 2: Details on SM background processes and their categorisation. Fragmentation, 
hadronisation and final state radiation were simulated with pythia 6.1 [T7j. photos [IB] 
was used to model final state radiation, and tauola [19j was used for tau decays. More 
details can be found in Ref. [20] . 



Category 


Process 


Software 




e + e~ -> Z/7* -»■ qq(7) 


KK 4.14 [21] 




Bhabha and e + e" -> Z/7* -> e + e"(7) 


BHWIDE 1.01 [22] 


2f 


e + e~ -»■ Z/7* -> ^+^(7) 


KK 4.14 [21] 




e + e _ -> Z/7* -> r + r _ (7) 


KK 4.14 gl] 




e + e~ — >■ Z — >■ ^(7) 


PYTHIA 6.1 [T71 




e+e" -> Z/7* ->■ W+W- 


KORALW 1.51 [23J 




e+e" -> ZZ 


PYTHIA 6.1 [T7] 


4f 


e + e~ — >• Z e + e~ 


PYTHIA 6.1 [17] 




e + e~ — )■ Z z/z/ 


PYTHIA 6.1 [T7] 




e + e" ->■ W ± e =F z/ 


PYTHIA 6.1 [171 


77 


77 -»■ 


phot02 [211 [25] 


77 -»■ qq 


phot02 [22125] 


n7 


e + e~ — >■ n7 


PYTHIA 6.1 [17] 



4 Event Selection 

A detailed description of the event selection criteria for the different Z decays, namely 
Z — > e + e~, Z — > fi + fi~ , and Z — >• i/p is presented below. The event topology is discussed 
first, followed by an explanation of how the visible decay products of the taus were treated. 
The Z reconstruction algorithm is then described and, finally, a detailed list of the cuts 
employed in the analysis is given. 

For the mass range considered, the Higgs boson is produced approximately at rest, 
and thus the decay h — > 2a — > 4r results in a pair of taus recoiling against another pair of 
taus. The JADE algorithm J2S1I27] was employed to cluster into jets all energy-flow objects 
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except for those identified as energetic, isolated photons, energy deposits in the lcal and 
SICAL, and the two hardest, oppositely-charged leptons in the case of the Z — > £ + £~ 
channels. Given that each jet is expected to arise from the on-shell decay a — > t + t~ , an 
effective way to target the signal topology is to use the jade algorithm with y cut chosen to 
merge proto-jets up to a mass of m- ]et = 15 GeV/ c 2 . This technique allows for an accurate 
clustering of the two tau pairs separately, thus rendering individual tau identification 
unnecessary and making track-multiplicity an effective discriminator against background. 

Because the taus from the same a decay are highly collimated, the identification of 
jets containing the decay products of two taus was based only on the track multiplicity 
of the jets, denoted nf ack , with the index i ordered in decreasing jet energy. Because 
the tau predominantly decays either to one charged particle ( "one-prong" decay) or three 
charged particles ("three-prong" decay), each jet is expected to contain two, four, or six 
tracks. To maximize the tracking efficiency, the jets were required to be well contained 
in the tracking volume. No distinction was made between leptonic and hadronic decays. 

The Z — > decay is often accompanied by additional photons from final state 

radiation, which can carry substantial momentum. An object is considered as an isolated 
photon if it is identified as a photon by the energy-flow algorithm, has E > 10 GeV, and 
contains less than 5% of the visible energy of the event in a cone of 10° around it. The 
photon was considered part of the candidate Z system when the invariant mass of the 
system was closer to the Z mass than the invariant mass of the lepton pair alone. 
This algorithm resulted in an increase of ~ 20% in the signal efficiency after the Z mass 
window cut, 80 < m z < 102 GeV/c 2 . 

For each of the channels below, a loose selection and final selection are presented. 
The loose selection isolates the broad characteristics of the signal events and allows for 
comparison of the data and simulated backgrounds. 

4.1 Z £+£- 

The loose selection consisted of the following requirements. An e + e~ or pair and 

the presence of two jets (or 3 jets with n| rack < 2 ) were required for consistency with 
the final state of the signal. The three-jet events are kept to recover signal efficiency for 
events with converted photon arising from final state radiation. Proper containment of 
the jet in the tracking volume was ensured by requiring | cos^jj < 0.9 and | cos#j 2 | < 0.9, 
where 9- h is the angle of the i th jet with respect to the beam axis . Additional lepton 
isolation was imposed by requiring that a cone of 10° around each lepton contained less 
than 5% of the visible energy of the event and cos^jj 1 " 1 < 0.95, where 9™ n is the minimum 
angle between each pairing of a jet and lepton. 

The final selection consisted of the following requirements and maintained an 
acceptable signal efficiency while rejecting most backgrounds. A mass window for the 
candidate Z between 80-102 GeV/c 2 was effective at removing two-fermion backgrounds. 
Due to the neutrinos from tau decays the signal was separated from fully hadronic final 
states by requiring $ > 20 GeV, where $ is the missing energy in the event. The 
expected jet configuration of the signal was enforced by requiring cos^^ < 0, where 9- n - i2 
is the angle between the two jets. Finally, the remaining backgrounds were suppressed 
by requiring n^2 Ck = 2 or 4, the dominant track multiplicities expected in the signal. 
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Figures [2] and [3] show the distribution of the reconstructed Z mass and missing energy 
for the Z — > e + e~ and Z — > fi + {i~ channels, respectively. The numbers of events passing 
loose and final selection in data and simulated background are shown in Table [3j 

4.2 Z vv 

All objects found in the event were clustered into jets as described above. The loose 
selection consisted of the following requirements. Missing energy greater than 30 GeV 
and missing mass, rfi, greater than 20 GeV/c 2 were used to reject dijet and other two- 
fermion backgrounds. In order to further reject the 77 background, events were required 
to have E vis > 0.05 E CM and | cos# me | < 0.9, where E vis is the visible energy and 9 me is the 
angle between the missing momentum vector and the beam axis. Events were required to 
have two well-contained jets with | cos^l < 0.85, dijet invariant mass m^ 2 > 10 GeV/ c 2 , 
dijet angular separation cos 6^,, < 0, and the highest energy jet was required to have 
E h > 25 GeV and < ack = 2 or 4. 

The final selection consisted of the following requirements. First, the requirement 
E- n + E- n + $ > E CM — 5 GeV was used to reject events with energy deposits in the 
forward regions of the detector. Consistency with Z — y vv was ensured by requiring 
$ > 60 GeV and rfi > 90 GeV/ c 2 . The distribution of aplanarity for the signal is 
strongly peaked near 0, while the remaining backgrounds have a longer tail. The tail 
of the aplanarity distribution for the signal extends further for larger m a and smaller 
m h ; larger m a leads to broader jets and lighter Higgs bosons can be produced with more 
momentum reducing the opening angle between the jets in the laboratory frame. Thus 
the requirement aplanarity< 0.05 was chosen to maintain an acceptable signal efficiency 
for mh = 86 GeV/ c 2 and m a = 10 GeV/ c 2 . Finally, the second jet was also required to 
have ri2 rack = 2 or 4. Figure H] shows the distribution of dijet invariant mass and missing 
mass for the Z — > vv channel. The numbers of events passing loose and final selection in 
data and simulated background are shown in Table |3J 

4.3 Signal Efficiency 

The h — > 2a — > 4r signal efficiency is shown in Fig. as a function of the Higgs boson 
mass with m a = 4-10 GeV/ c 2 for the three Z decay channels considered. The decrease 
in efficiency at higher m a values, seen in Fig. El is due to two effects. First, the invariant 
mass of the jet becomes larger and the fraction of three-jet events increases. Second, and 
more importantly, the events become more aplanar. 

5 Systematic Uncertainties 

Uncertainties and inaccuracies in the Monte Carlo simulation lead to systematic effects 
in the analysis. The impact of uncertainties in jet energy and direction, missing energy, 
and lepton identification and isolation were estimated in Ref. [2B] for similar final states. 
Compared to those, the present analyses do not use neural networks for event selection, do 
not use b tagging or tau tagging, and the simulated background samples are substantially 
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Table 3: Number of events passing loose and final selections in each channel, in data, 
simulated background, and simulated signal (mh = 100 GeV/c 2 , m a = 4 GeV/c 2 ). The 
numbers of events passing the final selection are categorised by track multiplicity. 



Channel 


Selection 


data 


total 


background category 


signal 




track y~, track \ 

vn j u 2 ) 




background 


Of 
LI 


41 


77 








Loose 


299 


o o o 

332 


183 


137 


12.31 


0.65 


2.27 


Z -> e+e" 


(2,2) 
(2,4) + (4,2) 






0.034 
0.055 


0.034 
0.014 


0.000 
0.005 


0.000 
0.037 


0.000 
0.000 


0.689 
0.610 




(4,4) 





0.031 


0.019 


0.013 


0.000 


0.000 


0.126 




Loose 


83 


74.50 


12.79 


60.64 


1.07 


0.00 


2.37 


Z — > n + n~ 


(2,2) 
(2,4) + (4,2) 






0.058 
0.005 


0.005 
0.000 


0.053 
0.005 


0.000 
0.000 


0.000 
0.000 


0.800 
0.676 




(4,4) 





0.006 


0.000 


0.006 


0.000 


0.000 


0.127 




Loose 


206 


200 


135 


47.97 


13.50 


3.74 


12.63 


Z — > vv 


(2,2) 
(2,4) + (4,2) 






1.312 
1.948 


0.663 
0.528 


0.408 
0.575 


0.240 
0.845 


0.000 
0.000 


5.097 
4.741 




(4,4) 


2 


2.569 


0.461 


0.820 


1.288 


0.000 


1.089 



larger. The systematic uncertainties related to the simulation of jet energies and directions 
were evaluated from the sample of hadronic events collected at the Z peak in 1998. Based 
on that sample, it was found that additional smearing of the Monte Carlo simulation was 
not necessary in the barrel region of the detector. 

For the Z — > £ + £~ channels, the total relative systematic uncertainties from lepton 
identification and isolation were found to be 0.6%, 2.6% and 7.5% for the signal, ZZ, and 
Zee backgrounds, respectively. The systematic uncertainties for WW, Wei/, qq, and other 
backgrounds were all smaller than 30%. Based on these estimates and the background 
composition, a 10% uncertainty is estimated for the background in the Z — > channels. 

The cuts used for the Z — > vv final state are sensitive to beam related backgrounds. 
The energy distribution of this background was measured with events recorded at random 
beam crossings. Additional energy depositions at angles below 12° were added randomly 
to all simulated events according to this energy distribution. The relative uncertainty 
in the total selection efficiency for the analyses presented in Ref. [28] was 5% for the 
signal and 10% for ZZ, and it is between 30% and 100% for the other background 
processes. Based on these estimates and the background composition, the uncertainty 
for the background in the Z — Y vv channel is estimated to be 30%. 

The agreement between the background estimate and the observed number of events 
in data with the loose selection is within the systematic and statistical uncertainty for 
all three channels. Given the low numbers of selected events, the final measurements are 
statistically limited. 
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6 Results 



No excess of events above the background was observed. Limits on the cross section 
times branching ratio with respect to the SM Higgsstrahlung production cross section, 

— > aa) x B(a — > t + t ) 2 were determined as follows. A joint 
probability density was constructed to describe the number of events in each of three jet- 
multiplicity pairings (indexed by m) for each of the three final states (indexed by /). The 
three jet-multiplicity pairings correspond to the different permutations of one-prong and 
three-prong tau decays in each of the jets, ignoring those with six tracks in an individual 
jet, leaving the three permutations (nf ack , n2 rack ) G {(2, 2), (2, 4) or (4, 2), (4,4)} = M. 
The event count N m j in each of these nine categories was modeled with a Poisson 
distribution about the sum of the uncertain background b m j and the expected signal 
s m j scaled by £ 2 . A normal distribution was used to model the relationship between 
the uncertain background, the Monte Carlo-based background estimate b^j, and its 
systematic uncertainty A/. This procedure leads to the following joint probability density 
for the event counts: 

P(N mJ \e,b mJ )= n n Pois(N mJ \es mJ + b mJ ) ■ N(b™ c f \b mJ ,A f ). (1) 

m£M f£{ee,fj,fj,,i/u} 

Confidence intervals were constructed by using a generalized version of the Feldman- 
Cousins technique |29j, which incorporates systematic uncertainties in a frequentist 
way [30j EH]- Figure shows the 95% confidence level upper-limit on £ 2 as a function 
of for m a = 10 GeV/c 2 . Figure Eb shows 95% confidence level contours of £ 2 in the 
(m h ,m a ) plane. Because the selection has no m h or m a dependence, the resulting upper 
limits are fully correlated. The observed number of events is consistent with a downward 
fluctuation of the background, which leads to stronger than expected limits on £ 2 . 

7 Conclusions 

A search for a Higgs boson produced via Higgsstrahlung at LEP2 energies has been 
performed, where h — > 2a — > 4r and Z — > e + e~,fi + fi~, vv. No evidence for an excess of 
events above background was observed, and a limit on the combined production cross 
section times branching ratio, £ 2 = g g^ e ( e + e ^zh) x ~~ ^ aa ) x ~~ ^ r+r ~) 2 is 
presented. For m h < 107 GeV/c 2 and 4 < m a < 10 GeV/c 2 , £ 2 > 1 is excluded at 
the 95% confidence level. This analysis covers a region of parameter space previously left 
unexplored, and further constrains models with non-standard Higgs decays, such as the 
NMSSM. 
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Figure 2: Distributions for the Z — > e + e~ channel after the loose selection for (a) the 
reconstructed Z invariant mass and (b) missing energy, where signal corresponds to 
mh = 100 GeV/c 2 , m a = 4 GeV/c 2 with £ 2 = 1 (see text). The same distributions 
are shown in (c) and (d) after the final selection, excluding any requirements on the 
variable shown. 
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Figure 3: Distributions for the Z — > fi + fi~ channel after the loose selection for (a) 
the reconstructed Z invariant mass and (b) missing energy, where signal corresponds 
to mh = 100 GeV/ c 2 , m a = 4 GeV/ c 2 with £ 2 = 1 (see text). The same distributions are 
shown in (c) and (d) after the final selection, excluding any requirements on the variable 
shown. 
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Figure 4: Distributions for the Z — > vv channel after the loose selection and requirement of 
1 < -n,2 iack < 7 for (a) dijet invariant mass and (b) missing mass, where signal corresponds 
to mh = 100 GeV/ c 2 , m a = 4 GeV/ c 2 with ^ 2 = 1 (see text). The same distributions are 
shown in (c) and (d) after the final selection, excluding any requirements on the variable 
shown. 
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Figure 5: Signal efficiency as a function of the Higgs boson mass for the three channels 
considered in this work, Z — > e + e~, and vv. The upper (lower) portion of the 

efficiency band corresponds to m a = 4 (10) GeV/c 2 . 




Figure 6: (a) Observed and expected 95% confidence level limit on £ 2 as a function of the 
Higgs boson mass for m a = 10 GeV/c 2 . (b) Contours of observed 95% confidence level 
limit on £ 2 in the (mh, m a ) plane. 
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